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We investigate the pairing symmetry of layered BiS2 compounds by assuming that electron- 
electron correlation is still important so that the pairing is rather short range. We find that the 
extended s-wave pairing symmetry always wins over d-wave when the pairing is confined between 
two short range sites up to next nearest neighbors. The pairing strength is peaked around the 
doping level x = 0.5, which is consistent with experimental observation. The extended s-wave 
pairing symmetry is very robust against spin-orbital coupling because it is mainly determined by 
the structure of Fermi surfaces. Moreover, the extended s-wave pairing can be distinguished from 
conventional s-wave pairing by measuring and comparing superconducting gaps of different Fermi 
surfaces. 
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Very recently, Mizuguchi et al discovered a new 
layered superconductor B14O4S3 with T c ~ 4.5 K[l| 
and other groups also reported similar materials, 
LaOi_ x F ;E BiS2[3 and NdOi_ x F :c BiS2i, 0. These ma- 
terials share some common features with both copper- 
based and iron-based high temperature superconductors. 
The electronic properties of these materials are deter- 
mined by the BiS2 layer [lj], similar to the Cu02 layer 
in cuprates or the Fe2An2(An=P, As, Se, Te) layer 
in iron-based superconductors. Electron doping to the 
BiS2 layer can induce superconductivity. For example, 
in LaOi_ x F a; BiS2, the partial replacement of O by F 
provides electron doping and induces superconductivity. 
Some properties in the superconducting (SC) state are 
also shown to be unconventional although it is too early 
to make a conclusive statement [5]. A number of stud- 
ies appear shortly after the discovery of BiS2 layered 
superconductor [6|48[ . 

The basic band structure of these materials has been 
calculated by first principle calculation. The conduc- 
tion band in the BiS2 layer is mainly attributed to the 
p x and p y orbits of Bi[8|]. Neglecting the interlayer 
coupling and states away from Fermi surface, Usui et 
a/.[8[ proposed a two-orbital model, in which only p x 
and p y orbitals are included, to describe the band struc- 
ture. A good Fermi-surface nesting and two quasi-one- 
dimensional bands have been found in LaOi_ x F ;E BiS2[l]. 
The importance of Fermi surface nesting has been em- 
phasized in Ref. 0. 

Because of the low T c and the extended p-orbitals 
in BiS2 layered materials, electron-phonon coupling was 
suggested to play the main role in the Cooper pairing [8, 
lllj . However, because electron-electron correlation gen- 
erally is more important in a low dimensional system, 
the correlation effect might play an important role in 
driving superconductivity even if the p-orbitals of Bi are 
much less localized compared with d-orbitals in cuprates 



and iron-based superconductors [8L_J10|. Experimentally 
it has also been proposed in Ref. [5| that the SC pairing 
is strong and exceeds the limit of the phonon mediated 
picture. The correlation effect, therefore, seems to be a 
good candidate responsible for the SC pairing in these 
materials. 

In the two-orbital model based on the p-orbitals of Bi, 
the quasi-one dimensional bands are generated through 
the p-orbitals of S[8]. If electron-electron correlation is 
important, the superexchange mechanism can naturally 
result in a next-nearest-neighbor (NNN) antiferromag- 
netic (AFM) exchange coupling. Although the existence 
of AFM fluctuation is still needed to be justified exper- 
imentally, theoretically it is legitimate to search for the 
consequence of possible correlation effect in these mate- 
rials. 

In this paper, we investigate the pairing symmetry 
of layered BiS2 compounds by assuming that electron- 
electron correlation is still important so that the pairing 
is rather short range. Under the assumption that the 
short range pairing stems from short AFM exchange cou- 
plings, we find that the extended s- wave pairing symme- 
try always wins over d-wave. Such result is very similar to 
the case in iron-based superconductors 12l-[l4| . We find 
that the pairing strength is peaked around the doping 
level x — 0.5, which is also consistent with experimental 
observation. The pairing symmetry is very robust against 
spin-orbital coupling because it is mainly determined by 
the structure of Fermi surfaces. The extended s-wave 
pairing can be tested by measuring the gap distributions 
at different Fermi surfaces. 

Theoretical Model. LaOBiS2 is an indirect semicon- 
ductor. Its band gap is about 0.82 eV, which is much 
larger than the superconducting gap, T c ~ 10 K. The 
simplest effective model is a two-orbital model con- 
structed by 6p x and 6p y orbitals of Bi. As shown in 
band structure of LaOo.5Fo.5BiS2 from first principle 
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calculation^, the four lowest conduction bands are split- 
ted slightly. The splitting is so small compared with their 
band width so that the low energy physics can be cap- 
tured by a two dimensional two-orbital model that ig- 
nores the coupling between BiS2 layers. We adopt the 
band structure proposed in Ref. HI, in which the nonin- 
teracting Hamiltonian reads Hq = a ^ a T(k)^k,cr, 



T(k) = 



e x (k) - fi e XY (k) 
e XY (k)* e Y (k)-fi 



(1) 



where ^ G — (c£ a x , a y) are creation operators for 
electrons with spin <j in the two orbitals Px,Py aligned 
along diagonal directions in Bi unit cell and 

ex(k) = 2t(cos(k x ) + cos(ky)) + 2t r cos(k x + k y ) 
-\-2t" cos(k x — k y ) 

+2t 2 i(cos(2k x + fc^) + cos(k x + 2fc y )) 
+2t / 21 (cos(2A: ;E — fc y ) + cos(— fc^ + 2k y )), 
ey (k) = 2£(cos(& :E ) + cos(/^)) + 2t' cos(/c ;E — fc^) 
+2t /; cos(^ + fc y ) 

+2t 2 i(cos(2k x — k y ) + cos(— fc^ + 2/c^)) 
-f 2t / 21 (cos(2£; ;E + k y ) + cos(k x + 2k y )), 
e X y(k) = 2t X y(cos(/c :c ) - cos(/c 2/ )) 

+2t 2 oxy(cos(2/c :c ) - cos(2& y )) 
+2t2ixy(cos(2/c ;E + fc y ) — cos(k x + 2/c y ) 
+ cos (2/^ — fc^) — cos(/c ;E — 2k y )). (2) 

The hopping parameters in Eq.[2]are given by t = —0.167, 
t' = 0.88, t" = -0.094, t xy = 0.107, t 21 = 0.069, 
t 21 = 0.014, t 2 oxv = -0.028, and t 2 ixY = 0.02 in units 
of eV. It is clear from these parameters that the material 
is quasi-one dimensional because the value of t' is almost 
an order of magnitude larger than other hopping param- 
eters. The model has a Lifshitz transition at electron 
doping x = 0.452 and x = 0.515 where the Fermi surface 
topologies are changed. 

If the electron-electron correlation effect can not be 
ignored, the hopping through the p-orbitals of S would 
naturally lead to an AFM exchange coupling between 
two NNN sites of Bi for the above two-orbital models. 
Taking the standard approach, in general, we can write 
the interacting Hamiltonian as 
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where S ? - 



\Y^oo' C i<TcP™ ,c w'& is tne local s P in 



notes a pair of NNN sites, a is orbital index and we 
use X, Y to label px and py orbitals, respectively. Ow- 
ing to the quasi-one dimensional property, the AFM 
NNN coupling can take two independent values such that 
Ji,x = J 2Y = an d J 2 x = ^2,y = J 2 , keeping C4 



symmetry unbroken, with J 2 expected to dominate over 
other exchange parameters. If we rewrite the interaction 
term in momentum space, 



k,k' 



N 



(J 2 , a cos(k x + k y ) cos(k f x + k' y ) 



(4) 



In general, the presence of electron-electron correlation 
can also significantly renormalize the bare band struc- 
ture. In the strong coupling limit, the band width renor- 
malization strongly depends on doping. However, here 
since the correlation is, at most, moderate, we assume 
that the renormalization does not vary significantly as 
the function of doping. Under such an assumption, we 
can use standard mean-field approximation to obtain the 
SC state and its favored pairing symmetry. This ap- 
proach has been used to study iron-based superconduc- 
tors where the correlation is most likely around interme- 
diate coupling strength. Under mean-field approxima- 
tion, the total Hamiltonian can be simplified as 

h mf = 4^(k)^ k + Yl ( e *(- k ) + e ^(- k ) - 



-—(At A x 
2J 2 V X 



AyAy + A X A X + Ay Ay), 



(5) 



where 



A(k) 



( e x (k) e XY (k) A xx (k)* \ 

e XY (k)* ey(k) Ayy(k)* 

A xx (k) -ejr(k) -e XY (k)* 

\ Ayy(k) -e XY (k) ey(k) J 



with 4>\. = (C^ t)X ,C^ t y,C_k4 > x,C'_k4 > y) , e a =x, Y 



(6) 



- /I, dk',t,a 

A aa (k) = 
A„ = 



(C-k'^aCk' 



and 



A' 



A a cos(k x + k y ) + A^ cos(k x — I 
2 . 

-^r^2,a ^2 ^,t,a cos (K + k y) 
k' 

■Jf J 2,a dk ''t,« cos (K " k 'y) 



(7) 
(8) 

0) 



Since we only consider spin singlet pairing, there are 
two possible pairing symmetries. In above mean field 
approximation, the s-wave pairing is taken if Ax = 
Ay while the <i-wave pairing is cho- 



A' Y and A^ 



operator, tt^q, is the local density operator, de- sen if Ax 



"Ay, A' x 



-Ay. The above equa- 



tions can be solved self-consistently with standard ap- 
proach. By diagonalizing A(k) via an unitary trans- 
formation, ?/1"(k)A(k)i7(k), we obtain four Bogoliubov 
quasi-particle eigenvalues E\ = —E3 and E 2 = —E^ 
which are given by 
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£ m=1 , 2 (k) = ^\[¥x + ?y + 2kxy| 2 + A xx + A YY )± 

±V^ " 4 + 4, " A^) 2 + ^xv\ 2 [(ex + ?y) 2 + + A YY ] - 4A*xA Fy (e| F + ef Y ). 

(10) 



The self-consistent gap equations are 

k , rn 

UZ Wtm (k)U m>m (k)F[E m (k)] (11) 
^ 2 

^x(y) = ^ ~J^^2,x(Y) cos (k x - k y ) 

k,m 

^3(4), m (k)^i(2), m (k)F[S m (k)] (12) 

where is Fermi-Dirac distribution function, F[E) = 
l/(e E / kBT + l). 

The above equations can be solved numerically. In 
Fig. [H we plot the free energy of s- and d-wave pairing 
states in J2 — J2/J2 plane with 0.51 electron doping. In 
all plotting region, the extended s-wave is favored. The 
reason that the extended s-wave is favored over the d- 
wave can be roughly understood by the overlap between 
the absolute value of the pairing gap form factor with 
the Fermi surface topologies as shown in Ref. 15. For 
simplicity, we consider the case that only J 2 exists. We 
first assume the s- wave and d-wave have the same pairing 
amplitude, then the overlap between SC gap and Fermi 
surface is defined as the sum of the absolute values of 
the SC gap form factor over the whole Fermi surface. 
The pairing channel is the one that has largest overlap. 
Fig. [2] illustrates the SC gap distribution over the entire 
Fermi surface with Ax = 0.2 at doping x = 0.51. The 
gap distributions of s- and d-wave pairing states are the 
same except the region near the nodes of d-wave state. 
It's obvious that the overlap of s- wave pairing state is 
larger than that of d-wave pairing state. Our calculation 
also shows that the s- wave pairing does not co-exist with 
(i-wave pairing. As shown in Fig. [TJ the mixed s-wave 
and d-wave state has higher free energy than the pure 
s- wave state. 

The extended s-wave pairing is very robust against the 
variation of parameters. This is because the quasi-one 
dimensional nature in the band structure results in the 
dominant pairing channel which is mainly from the J2 
term. For example, even we increase J' 2l the amplitude 
of pairing induced by J2 is still much larger than that by 
J2, as shown Fig. [3] (a). Ax is five times of A' x when J2 = 
0.6, J 2 = J 2 and ten times when J 2 = 1, J 2 — h- There 
is also a critical value of J 2- SC state is only favored when 
J2 is larger than the critical value. As shown in Fig. [3] (a), 
while J 2 < 0.2, superconductivity disappears at x = 0.51. 



Furthermore, the optimal doping to achieve the highest 
T c is around 0.5. Therefore, given J2, we will obtain 
a superconducting dome with optimal T c at x ~ 0.5 as 
shown in Fig. [3] (b). These results are consistent with 
the experimental results. Experimentally, the optimal 
doping of LaOi_ x F ;r BiS2 is suggested to be ~ 0.5 per Bi 
atom[2] and - 0.3 per Bi atom for NdOi_ :c F :c BiS2[3. In 
the later case, the actual value may be higher because of 
the existence of impurity phases [H, 

The robustness of the extended s-wave is preserved 
even if we include other short range AFM couplings such 
as nearest-neighbor (NN) ones or spin-orbital coupling 
(SOC). In the mean-field level, the NN AFM coupling 
favors d x 2_ y 2 pairing symmetry. As shown in Fig. 31 as 
soon as J 2 is larger than the NN coupling Ji, the ex- 
tended s- wave dominates. Furthermore, the extended 
5-wave is also very robust against spin-orbital coupling. 
Because of the absence of p z orbital contribution near 
Fermi surfaces, the onsite SOC Hamiltonian is given by 
Hso = -iA/2E kff T(Ck.Crk CT - C Yk<7 C Xka ). By fit- 
ting the first principle band structure with SOC with the 
two-band model, we get A ~ 1.1 eV. However, even if this 
SOC appears to be strong, it only lifts the degeneracy 
at X and T points and makes no significant adjustment 
on Fermi surfaces. We find that the extended s- wave is 
hardly affected. An only visible effect is that the SOC 
slightly increase the optimal doping level when we repeat 
the above mean-field calculation in the presence of SOC. 

In order to differentiate the extended s-wave from the 
(i-wave or a conventional s-wave, we plot gap values 
around entire Fermi surfaces at different doping levels 
that are characterized by different topologies as shown in 
Figs. [2j [5]and[6l For the extended s- wave, while the SC 
gap is developed throughout the Fermi surfaces, the gap 
has significant variations. The variations are different in 
three different Fermi surface topologies. In particular, 
at doping levels 0.51 and 0.69, the gap value around the 
corners of the square-shape Fermi surface at the center 
of Brillouin zone is significantly smaller than the values 
on the other Fermi surfaces. This feature can allow us 
to distinguish the extended s-wave from the d-wave and 
a conventional s- wave. In the d-wave, as plotting in the 
same corresponding figures, there are gapless nodes at the 
above corners. In a conventional s-wave, the gap values 
should be more or less isotropic around Fermi surfaces. 
Experimental measurements by the standard scanning 
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FIG. 1: The free energy of s- and d-wave pairing states as a 
function of J 2 , J 2 at 0.51 electron doping. 
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(b) (d) 



FIG. 2: The SC gap distribution over the whole Fermi surface 
with Ax = 0.2 at 0.51 electron doping: (a) for s-wave pairing 
state and (c) for d-wave pairing state. The line width of 
Fermi surface is proportional to [ ] ~ 1 ? where e(k) is the 
eigenvalue of T(k) and d *j^ represents the derivative of e(k) 
along normal direction of Fermi surface. The density of states 
(DOS) in SC ground state: (b) for s-wave pairing state and 
(d) for d-wave pairing state. 



tunneling microscope (STM) or angle-resolved photoe- 
mission spectroscopy (ARPES) can provide a verdict. 

Summary. If the pairing in BiS2-based superconduc- 
tor is induced by the moderate electron-electron correla- 
tion, we show that similar to iron-based superconductors, 
the pairing symmetry is an extended s- wave, which is 
robust against doping and other parameter changes be- 
cause of the quasi-one dimensional electronic structure 
of the materials. As the extended s-wave state is fully 
gapped with a significant gap variations on different parts 
of Fermi surfaces, our prediction can be directly justified 
or falsified by future experimental measurements. 
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2012CV821400), as well as national science foundation 
of China ( Grant No. NSFC-1190024). 
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FIG. 3: (a) The amplitude of pairing gap induced by J2, 
Ax, and that by J 2 , A f x , as a function of J = J 2 = J 2 
at 0.51 electron doping, (b) Ax as a function of doping at 
J 2 0.6, J 2 = 
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FIG. 4: The amplitudes of different pairing gap components 
induced by Ji,J 2 with J 2 = at electron doping= 0.51. 
snnn is the s wave resulting from NNN AFM coupling, 
s x 2 +y 2 and d x 2_ y 2 from NN AFM coupling. 
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FIG. 5: The SC gap distribution over the whole Fermi surface 
with Ax = 0.2 at 0.25 electron doping and its corresponding 
DOS: (a),(b) for s-wave; (c),(d) for d-wave. 
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(b) (d) 

FIG. 6: The SC gap distribution over the whole Fermi surface 
with Ax = 0.2 at 0.69 electron doping and its corresponding 
DOS: (a),(b) for s-wave; (c),(d) for d-wave. 
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